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ABSTRACT 



In this paper we describe the first data release of the UltraVISTA near-infrared imaging survey of the COSMOS field. We summarise 
the key goals and design of the survey and provide a detailed description of our data reduction techniques. We provide stacked, sky- 
subtracted images in YJHK S and narrow-band filters constructed from data collected during the first year of UltraVISTA observations. 
Our stacked images reach 5a AB depths in an aperture of 2" diameter of ~ 25 in Y and ~ 24 in JHK S bands and all have sub-arcsecond 
seeing. To this 5<x limit, our A" s catalogue contains 216,268 sources. We carry out a series of quality assessment tests on our images 
and catalogues, comparing our stacks with existing catalogues. The ltr astrometric RMS in both directions for stars selected with 
17.0 < Zf s (AB) < 19.5 is ~ 0.08" in comparison to the publicly-available COSMOS ACS catalogues. Our images are resampled to 
the same pixel scale and tangent point as the publicly available COSMOS data and so may be easily used to generate multi-colour 
catalogues using this data. All images and catalogues presented in this paper are publicly available through ESO's "phase 3" archiving 
and distribution system and from the UltraVISTA web site. 



Key words, observations: 
structure of Universe 



galaxies: general - galaxies: high-redshift - astronomical data bases: surveys - cosmology: large-scale 



(N 
> 

OO 

in 

\o 

O 

> 

X 

C3 



1. Introduction 

The vital role of near-infrared (A ^ 1 - 2.5 yum) imaging surveys 
for advancing our understanding of galaxy evolution has long 



been recognised (|Cowie et al.||1990{ |Glazebrook et al. 



While optical surveys utilising large-format charge-coup 



1991 



ed de- 



vice (CCD) detectors were the first to enable the discovery of 



substantial samples of normal galaxies at redshifts z > 2 ( Steidel 



et al. 1996 1 Madau et al. 1996 ), it was already known that at least 
some galaxies at high redshift were either too old or too dust- 
obscured to be easily detected by rest-frame near-ultraviolet se- 
lection (Dunlo p" et al.|1996 Dey et al.|1996) >. In addition, even 
for apparently young UV-luminous galaxies, the value of us- 
ing near-infrared observations to sample the rest-frame optical 
light, more representative of the evolved mass-dominant stel- 
lar population, was understood and indeed demonstrated before 
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* Based on data products from observations made with ESO 
Telescopes at the La Silla Paranal Observatory under ESO programme 
ID 179.A-2005 and on data products produced by TERAPIX and the 
Cambridge Astronomy Survey Unit on behalf of the UltraVISTA con- 
sortium. 

** Catalogs are available at the CDS via anonymous ftp to cdsarc.u- 
strasbg.fr (130.79.128.5) or from http://cdsarc.u-strasbg.fr/ 
|viz-bin/qcat?J/A+A/S44/A156p 



the advent of multi-pixel near-infrared imagers ( Lilly & Longair 
[T984| . 

However, near-infrared surveys are a challenging proposi- 
tion for several reasons. Firstly, at near-infrared bandpasses the 
sky background is extremely bright; in the K s band, in AB mag- 
nitudes, it is typically 15 mag/arcsec 2 , which means that many 
short exposures must be combined in order to avoid detector sat- 
uration on the sky, greatly increasing overheads. Secondly, the 
sky background is time-variable, many magnitudes brighter than 
the faint astronomical sources of interest, and so must be care- 
fully subtracted from each image before scientific exploitation 
can take place. Lastly, conventional silicon CCDs are very in- 
efficient at near-infrared wavelengths, and a different detector 
technology must be employed which is an order of magnitude 
more expensive. In terms of sky footprint, near-infrared detec- 
tors have generally lagged behind optical detectors by approxi- 
mately a decade. 

Nevertheless, these challenges have been progressively over- 
come, and, following the pioneering work described above with 
early near-infrared arrays such as IRCAM on the UK Infrared 
Telescope (UKIRT) ( |McLean et alj]T986l >, the full potential of 
near-infrared surveys to clarify our view of galaxy evolution at 
Z - 1 - 3 began to be realised with the advent of larger for- 
mat infrared array cameras such as ISAAC on ESO's Very Large 
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Telescope (V LT) ( |Cimatti et al.|2002| |Labbe et al.|2003| |Franx 



et al. 2003 1. Meanwhile, the importance of near-infrared sur- 



veys for revealing dust-obscured star-forming galaxies was fur- 
ther enhanced by the discovery of significant numbers of dusty- 
enshrouded high-redshift star-forming galaxies at sub-mm wave- 
lengths ( |Hughes et al.|1998|[Scott et al.|2002] i. Around the same 
time the unique power of the deepest near-infrared imaging to 
conduct rest-frame ultraviolet surveys for galaxies at z > 6.5 
was first demonstrated using the NICMOS camera on the Hubble 
Space Telescope (HST; Bouwens et al.|[2004 Thompson et al. 
[2005] l 

Despite these impressive advances, the field-of-view offered 
by near-infrared cameras such as IRCAM, NICMOS and ISAAC 
was very small (a few arcmin 2 ), and it is only in the last half- 
decade or so that the introduction of genuinely large-format 
near-infrared array cameras has enabled efficient, deep near- 
infrared imaging of degree-scale areas of sky, allowing stud- 
ies of more representative volumes of the high-redshift universe 
(i.e. ^ 100 x 100 comoving Mpc). First WFCAM on UKIRT 
( |Casali etaT 2007 ), then WIRCam on the Canada France Hawaii 
Telescope (CF HT) ( |Puget et al.|2"004| , and NEWFIRM at NOAO 
( |Probs t 2004) have heralded a new era of major coordinated 
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near-infrared survey programmes (e.g. UKIDSS, Lawrence et al. 
J; NEWFIRM Medium-Band Survey, |van Dokkum et al. 

I); 



(2009)); WIRD S and associated near-infrared follow-up of the 
COSMOS field ( |Bielby et al.|20TT]|McCracken et al.|2010] >. This 
has led to a number of breakthroughs in extra-galactic astron- 
omy, including, for example, the study of the bright end of the 
galaxy luminosity function from z — out to z - 6 ( |McLure 
et al.|2009{|Cirasuolo et al.|2010|l, and the d iscovery of the most 
distant known quasar ( |Mortlock et al.|201 l| i. 

In addition, deep, wide-field near-infrared photometry cou- 
pled with high-quality optical surveys has enabled spectral en- 
ergy distribution (SED) fitting techniques to be pushed beyond 
z ~ 1.5. Near-infrared data play a key role in minimising the 
catastrophic failure rates in photometric redshift estimates and 
provides robust rest-frame visible flux determinations at z ~ 2 
( Ilbert et al.||2009[ ), enabling measurements of the evolution of 
the mass buildup in stars over a large fraction of the age of the 
Universe ([Drory et al.|[2003| |Arnouts et aT||200"7t |Ilbert et al.| 
20T0l|Caputietal.|2011| i. 



These efforts have now culminated in VISTA (Emers on &| 
|Sutherland]|20 1 0[ ) the first 4-m class telescope specifically de- 
signed to conduct wide-area near-infrared surveys and equipped 
with a large-format array camera, "VIRCAM" (Dalton 2006). 
Thanks to its large mosaic of 16 detectors, VIRCAM is currently 
the most efficient wide-field near-infrared survey camera in the 
world (around four times more efficient than WIRCam, and three 
times as efficient as WFCAM). It also has the benefit of being 
mounted on a telescope for which virtually all observing time is 
available for surveys, and for which observations are efficiently 
programmed in queue-scheduled mode. Inspired by the success 
of UKIDSS, ESO has implemented a coordinated multi-tier pub- 
lic survey programme with VISTA. The UltraVISTA survey pre- 
sented here is the deepest component of the VISTA survey "wed- 
ding cake". 

Covering an area of 1.5 deg 2 , UltraVISTA is significantly 
larger than the only comparably-deep near-infrared survey con 



the planned exploitation of UltraVISTA for the discovery of the 
most luminous galaxies at z — 7, e.g., Bowler et al. (2012). 

To maximise the leverage and legacy value of these new 
deep near-infrared data, the UltraVISTA survey is centred on the 
COSMOS field, the location of the largest ever ACS optical mo- 
saic obtained with HST (Sc oville et al. 2007; Koeke moer et al. 
2007) and an ever growing heritage of deep ground-based and 
space-based multi-frequency imaging and spectroscopy The 
first-year data set described in this paper is already deeper than 
all existing COSMOS NIR data ( |McCracken et al.|2010|[Bielby 
et al. 201 1\ in all bands by between one and two magnitudes and 
also contains for the first time deep Y— band imaging. 

To most efficiently exploit VISTA for the discovery and 
study of UV-selected galaxies at the highest redshifts (z — 
6.5 - 9) and in the investigation of the growth of galaxies 
through the crucial redshift range 1 < z < 3 when cos- 
mic star- formation density peaks (Hopkins & Beacom 2006), 
the UltraVISTA survey comprises three separate components: 
a wide, deep Y, J, H, K s survey (a contiguous field covering a 
1.5 deg 2 ); an ultra-deep Y,J,H,K S survey (consisting of deeper 
strips covering a 0.7 deg 2 ), and an ultra-deep narrow-band (A = 
1.1 8 /mi) survey targeting emission-line galaxies at a range of 
redshifts, e.g. Ha at z = 0.8, [OIII] -emitters at z = 1.4, [Oil] 
emitters at z — 2.2, and ultimately Lya emitters at z = 8.8. To 
accomplish these goals, UltraVISTA has been allocated 1800hrs 
of execution time. 

It is important to stress that while the advent of Wide Field 
Camera 3 (WFC3) on HST in 2009 has enabled extremely deep 
near-infrared imaging (up to A a 1.6 /urn) which has revolu- 
tionised the study of galaxies at z — 7 - 8, (|Bouwens et aL 
[20T0| [McLure et al pOlOl |Oesch et aI1[2009l |Fmkelstein et al. 
2010; Bunker et al. 2010l the very small field-of-view offered 
by WFC3/IR coupled with its inability to observe in the K- 
band means that deep ground-based surveys such as UltraVISTA 
remain of crucial importance. In particular, the largest cur- 
rent (or indeed planned) WFC3/IR extragalactic survey is the 
Cosmic Assembly Near-infrared Deep Extragalactic Legacy 



Survey (CANDELS; Grogin et al. (2011 1; Koekemoer et al. 



ducted to dat e (the UKIDSS Ultra Deep Survey(UDS; |Almaini 



(2011 1), but even this 900-orbit 3-year HST Treasury Program 
will only cover 800 arcmin 2 . Thus UltraVISTA is an excellent 
complement to CANDELS, and indeed CANDELS has recently 
completed deep J, //-band WFC3/IR imaging of a a 200 arcmin 2 
region within the 1 .5 deg 2 UltraVISTA imaging described here 
(i.e. covering only - 4% of UltraVISTA). 

In this paper we present a detailed description of the data re- 
duction methods and properties of the five near-infrared stacks 
created from the first season of UltraVISTA operations. Already, 
with only these first images, the UltraVISTA survey has the 
largest etendue of any near-infrared survey 

All magnitudes in this paper, unless otherwise noted, are 
given in the AB system. Data products described here are avail- 
able from ESCQ the UltraVISTA websi tq^and CESAM[3 

2. Observations and data reductions 

2.1. Observations 

The images described here were taken between 5th December 
2009 and the 19th of April 2010 with the VIRCAM instrument 



et al.| |2007| l), and will ultimately go significantly deeper. 1 I http : //cosmos . astro . caltech . edu 



VIRCAM also offers two significant advantages over WFCAM 
(and indeed WIRCam or NEWFIRM) in that its Raytheon de- 
tectors are much more sensitive in F-band, and are essentially 
free from the electronic cross-talk. These are crucial benefits in 



2 http : //www . eso . org/sci/observing/phase3/data_ 

releases.html 

i http : //www . ultravista . org/ 

4 http : //ce sam . oamp . fr/ultr avi st a/index . php 
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Red filled: UltroVISTA Ultra-Deep stripes 

Blue hatched: UltroVISTA Deep stripes 

Blue outline: UltroVISTA contiguous region (depth = Deep or Ultra-Deep) 

Green outline: HST/ACS 
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Fig. 1. Schematic layout of UltraVISTA observations, showing 
deep and ultra-deep regions (hatched and filled regions respec- 
tively). The data described in this paper correspond to a uniform 
coverage in YJHKs of the contiguous region and to NB1 18 ob- 
servations of the ultra-deep stripes. 



on the VISTA telescope at Paranal as part of the UltraVISTA 
survey program. VIRCAM is a wide-field near-infrared camera 
consisting of 16 2048 x 2048 Raytheon VIRGO HgCdTe arrays 
arranged in a sparse-filled array with gaps between each array of 
0.90 & 0.425 of a detector in X and Y respectively ( Emer son & 
Sutherland|2010| . The mean pixel scale is 0.34"pixel 1 (Dalton 
2006) . 



The sky coverage of the 16 non-contiguous detectors is 
called a "pawprint". A contiguous region of size 1.5° x 1.23° can 
be covered by means of six pawprints suitably spaced in right 
ascension and declination with random 60" jitter offsets in both 
directions (two w 0.1° bands at the top and bottom of the field 
receive half the exposure time). 

Specifically, three pawprints with identical RA and with 
Dec differing by 5.5' = 47.5% of a detector height make up 
a set of four stripes (corresponding to the ultra-deep stripes in 
UltraVISTA), and another three pawprints shifted by 95% of a 
detector width in RA make up another set of stripes, which to- 
gether form a contiguous region where most pixels in the result- 
ing stack are covered by two of the six pawprints. 

Fig. [T] illustrates the layout of UltraVISTA observations 
showing the deep survey, which will cover the full survey area, 
and the ultra-deep part, which covers half of this area in a series 
of ultra-deep stripes. The first season of UltraVISTA data de- 
scribed in this paper comprises six contiguous pawprints in four 
broad-band filters covering the deep survey area, each with equal 
exposure times, and narrow band observations on the ultra-deep 
stripes; subsequent observing seasons are expected to concen- 
trate exclusively on the ultra-deep stripes. 

The observations, carried out in service mode, are specified 
by observation blocks (OBs). The characteristics of the OBs 
used in UltraVISTA season one are listed in Table [T] Most of 
the season one OBs comprise images jittered around the centre 
of a single pawprint position, with the jitters being drawn from 
a random, uniform distribution over a box of side length 120" 
(random jitters are necessary because of persistence effects in 
VIRCAM and are also essential to derive a good sky frame). 



The exception to this was the "NB118 three paws" OBs 
(Table [T] which comprised images jittered around the centres of 
the three pawprints forming the ultra-deep stripes. For OBs con- 
taining more than a single pawprint per OB, the nesting (Table[T]i 
is important, and we did not use the optimal value. These OBs 
had a nesting of "FJPME" such that F (filter) is the outermost 
loop, and E (expose) is the innermost loop. The important as- 
pect here is that the three pawprints (P) (spaced exactly by 5.5' 
in Dec) are completed before a random jitter (J) is applied. This 
means that the faint persistent images (i.e. fake sources that are 
memories of a bright star at that x,y position on the detector in 
the one or two previous exposures) will be present in the stack at 
positions located 5.5' (and 11') away from bright stars in DEC. 
We deal with this by masking the persistent images in the in- 
dividual NB118 images (see Milvang-Jensen et al., in prep, for 
details of the procedure). For the other UltraVISTA OBs, the 
faint persistent images are fully removed by the sigma clipping 
used in producing the stacks, thanks to the random jitters applied 
between each single exposure. The first season of observations 
described here comprise around 200 OBs in total. The average 
efficiency (calculated as the total exposure time divided by total 
execution time these OBs) was 77%. 

In light of our experience gained in the season one observa- 
tions described here, from season 2 onwards we modified some 
of the OBs. For Y, we changed the DIT to 60 sec (with NDIT 
= 2), since 30 sec was unnecessarily short; for H, we changed 
the DIT to 10 sec (with NDIT = 6), for the same reason. For 
NB1 18, we changed the DIT to 120 sec (with NDIT = 1), since 
300 sec was unnecessarily long. We also changed our observa- 
tion strategy to jitters centered around a single pawprint per OB, 
and changed the total exposure time per OB to 1 hour (corre- 
sponding to 30 jittered exposures in an OB). 



2.2. Image selection and grading 

VIRCAM images are transferred to the Cambridge Astronomy 
survey Unit |(CASU f1for pre-preprocessing and removal of the 
instrumental signature. This includes dark subtraction, correc- 
tion for rest anomaly, flat-fielding, initial sky-subtraction, de- 
striping, non-linearity corrections and gain normalisation( |]jwi"n| 
2004). CASU subsequently provides these pre-processed images 
for each survey, as well as stacks of images from a single OB and 
pawprint, comprising typically 30 or 60 images. 

For UltraVISTA we start from the individual pre-processed 
images, rather than the stacked OB blocks, for a number of rea- 
sons: firstly, the OB blocks are combined at CASU at the native 
pixel scale of the instrument, which means that in good seeing 
conditions (median FWHM ~ 0.6") VIRCAM data is under- 
sampled. For this reason it is preferable to re-sample these data 
at a finer pixel scale; secondly, one of the principal scientific 
aims of the UltraVISTA project is to make measurements of dis- 
tant (z > 6) and faint (K s ^ 24) galaxies. To do this requires 
extremely accurate removal of the sky background for each in- 
dividual image; in the version of the CASU pipeline we used, a 
single sky background was used for all images coming from a 
given OB, and objects were not masked using the deepest pos- 
sible mask. Given that the sky background is known to vary on 
shorter timescales, this process may lead to a systematic mag- 
nitude offset at faint magnitudes near bright sources. For these 
reasons we use an iterative sky-background removal technique 



5 |http://cas u . ast . cam . ac . uk/ surveys-pro j ects /vista/ 
technical/data-processing 
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Table 1. Characteristics of the OBs used in UltraVISTA season 1 



OB description 


Filter 


JJll Lsj 


MmT 
JNtJii 


Jitter parameters 
Pattern Amplitude ["] 




Nesting N pawprillts 


Total exp. [s] 




Y 


Y 


30 


4 


Random 


60 


30 


FPJME 1 


3600 


36 


J 


J 


30 


4 


Random 


60 


30 


FPJME 1 


3600 


37 


H 


H 


6 


10 


Random 


60 


60 


FPJME 1 


3600 


36 


K, long OB 


A', 


10 


6 


Random 


60 


60 


FPJME 1 


3600 


18 


K s short OB 


K, 


10 


6 


Random 


60 


30 


FPJME 1 


1800 


27 


NB 118 single paw 


NB118 


300 


1 


Random 


61 


11 


FPJME 1 


3300 


6 


NB 1 1 8 three paws 


NB118 


280 


1 


Random 


61 


4 


FJPME 3 


3360 


4 



Notes - The "Amplitude" column gives the Maximum Jitter Amplitude, where a value of 60" corresponds to jitter positions being drawn from 
a random, uniform distribution over a box of side length 120", centered on the nominal centre coordinates of the given pawprint. The "Nesting" 
column indicates the order in which different operations are done, see text. The "N pawpl . ints " column gives the number of pawprints done by the given 
OB. To cover the contiguous UltraVISTA field in an approximately uniform manner, a set of 6 OBs of type N pawpr j n t s = 1 are needed, each centered 
on the pawprint in question. The "Total exp." column gives the total exposure time contained in the OB; this number is DITxNDITxNj it xN pawprints . 
The "Nqbs" column gives the number of OBs of the given type that are associated with the data covered by this paper. 




seeing seeing seeing ellipticity ellipticity ellipticity 



Fig. 2. Seeing (left) and ellipticity (right) distributions for all UltraVISTA images considered. The arrow represents the median of 
each distribution for images classified as A or AB. Note that the distributions for each quality class have been re-normalised, and 
the vertical axis has been rescaled. A, AB, C represent the quality classifications described in the text. 



starting from the pre-processed images and also resample all 
data to a pixel scale of 0.15"pixel 

The images in this release were taken between 5th December 
2009 and the 19th of April 2010. This does not consist of the 
complete number of images taken for the UltraVISTA program 
in the 2009-2010 observing season; subsequently, around 10% 
additional images in H and Ks were made available by CASU 
using a different pipeline processing, after we had already graded 
the first batch of images. In order to maintain as a homogenous 
as possible data set, we restrict in this release ourselves to this 
initial batch. However, had we included these data, the average 
exposure time per pixel would have been 4800 s and 5400 s 
higher in H and K s respectively, i.e. only around 10% larger. 
The images considered here were all processed with v0.8 of the 
CASU pre-processing pipeline and in total, we consider 7031 
individual images (each of which is a single multi-extension fits 
image containing 16 image extensions one for each of the VISTA 
detectors). 

Since UltraVISTA represents the first significant amount of 
data from VIRCAM processed at TERAPIX, we wished to visu- 
ally inspect all images to identify any problems which had been 



potentially overlooked by the automatic pipelines. Therefore, all 
images were inspected and graded in the YOUP^] environment. 
Images were assigned a grade of A , B (usable for science), C, 
D (rejected). The left and right panels of Fig. [2] shows the see- 
ing FWHM (measured assuming a Gaussian core), ellipticity and 
grading distributions for all images. Based on these distributions 
we decided to keep all images which have stellar FWHM < 1 .0" 
and ellipticity < 0. 1 and which were classified as either A or B 
based on visual inspection. The visual inspection process in gen- 
eral finds images which have bad PSFs or other optical defects 
which would have not been found by a typical seeing or ellip- 



http : //youpi . terapix . fr/ 
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ticity cuQ In total we reject 426 images or around 6% of the 
total. 

We do not use the confidence maps provided by CASU, but 
create our own weight maps from the supplied flat-fields and bad 
pixel maps using the weightwatcher tool ( |Marmo & Bertin 
2008). For our NB118 images which were taken at a fixed set 
of jitter patterns and thus suffer from image persistence effects, 
we mask the persistent images using the procedure described in 
Milvang-Jensen et al. (in prep.). 

2.3. Two-step sky subtraction 



To derive our sky-subtracted images, we use a set of tools devel- 
oped at TERAPIX which run under the distributed processing 
environment "condor' fj (Th ese processing steps are described 
fully in Bielby et al. ( |2011 1). Sky-subtraction is a two-step it- 
erative process. To summarise, we start by adding back the sky 
background frames subtracted by CASU (which are supplied as 
part of the original data release.) Based on the first-pass stack 
(computed using the CASU sky-subtracted images) and astro- 
metric solutions, we compute object masks for each individual 
image. Next, we use these object masks (appropriately resam- 
pled based on an initial astrometric solution) to effectively re- 
move objects computed from a running sky for each individ- 
ual image, based on a median of images taken during a 20- 
minute sliding window. After the subtraction of the running sky, 
we re-"destripe" the images and remove large-scale background 
gradients using sextractor (Bertin & Arnouts 1996). In gen- 
eral, computing sky backgrounds for each of the 7000 images is 
highly processor intensive; for each image, it takes around 15-20 
minutes on a standard TERAPIX computing node. 



2.4. Astrometric and photometric solutions 

After sky-subtraction, weight maps and catalogues are computed 
once more for each image using QualityFITS. Saturated ob- 
jects, based on an examination of the distribution of objects in 
the peak surface brightness / magnitude plane, are flagged in 
these catalogues, and the weight maps are used to flag cosmic 
rays. Bad pixels are also flagged. Next, these catalogues are used 
to compute the final astrometric and photometric solutions which 
will be used to combine and scale the images. Astrometric so- 
lutions are computed independently from each filter using the 
scamp tool ( |Bertin| 2006), but use a common astrometric refer- 
ence catalogue drawn from the COSMOS ;-band CFHT data 
(the same reference catalogue used in Capak et al. ( 2007 1 and 



McCracken et al. ( 2010) >). We use a third-order polynomial solu- 



tion in x and y detector co-ordinates (note that unlike the CASU 
reductions, we do not assume a radially symmetric astrometric 
solution). In order to derive a more robust astrometric solution, 
we use a precomputed " . ahead" file for all images which spec- 
ifies the relative positions and orientations of each of the six- 
teen detectors. In addition, we require that all the detectors share 
a common tangent point (focal plane mode "SAME_CRVAL" in 



7 Some of these bad PSFs were caused in part by software errors in 
early versions of ESO's Survey Area Definition tool: all season 1 OBs 
had pointing centres such that when a jitter jump went too far in one 
direction, the guide star fell outside the guide CCD; guiding was not 
active for the remaining images of that OB. This was fixed in season 
two observations by moving the pawprint centers. These tracking er- 
rors produce double-lobed PSFs in some images; each of the individual 
PSFs are smaller than the requirement and so pass our cut. 

8 http://www.cs.wisc.edu/condor/ 



scamp). These steps ensure that we can reliably match our refer- 
ence astrometric catalogues for many thousands of input images 
(note that we do not use the higher-order terms of the initial as- 
trometric solution provided by CASU^J. Thanks to our densely 
sampled astrometric reference catalogue the internal sigma of 
our astrometric solution is ~ 0.08", 0.09" in directions North- 
South, East-West directions respectively. 

Compared to our reference catalogue, in the same directions, 
we find standard deviations of ~ 0.09", 0.10". Similar values are 
found in all filters. Given that the native pixel scale of VIRCAM 
is 0.34"pixeF 1 , our astrometric solution is more than sufficient 
to provide a precise and reliable image coaddition (in fact, our 
astrometric accuracy is probably limited by undersampling in the 
VIRCAM images). 

Our initial magnitude zero points for each individual im- 
age are based on those supplied by CASU for their .st stacks 
(which comprise a stack of several individual images), which is 
based on their calibration of the VISTA photometric system's 
zero points. To account for possible photometric variations be- 
tween the images in each .if stack we calculate a rescaling factor 
for each using scamp based on overlapping paw-prints. Note that 
the same rescaling factors are applied to all detectors: we assume 
that the relative scaling factors between chips does not change 
(the CASU processing pipeline equalises the gain between all 
detectors at the flat-fielding stage, and should remain constant). 
To create our final stacks in the AB magnitude system ( |Oke| 
[1974 1 we simply apply the appropriate flux scaling to convert 
the supplied Vega magnitudes to AB, based on the VISTA tele- 
scope detector, filter and atmosphere combination. The conver- 
sion factor C from AB to Vega we use are as follows, in the sense 
mag AB = mag vega + C where C = 0.61,0.90, 1.38, 1.84,0.86 for 
Y, J, H, K s and NB 1 1 8 filters respectively. 

Note that CASU produces "flat" images which have constant 
flux per pixel for a uniform illumination; this is taken into ac- 
count in the resampling stage. 



2.5. Coadded images 

In the last processing step, the images and weight maps are coad- 
ded using a modified version of the swarp software ( |Bertin et al.| 
2002 ) which permits a combination of images based on a clipped 
sigma estimator; we use a clipping threshold of 2.8<x. Before 
stacking, a small number of images which have large photo- 
metric extinction or bad astrometric solutions are also rejected. 
For the final stacks, in the five bands, 6520 images were used. 
Since the size of VIRCAM pixels varies radially as a function of 
distance from the centre of the mosaic, this must be accounted 
for during image co-addition. Bad regions on individual detec- 
tors (such as half of detector 16, whose pixels suffer from time 
variable quantum efficiency, most notable at shorter wavelengths 
where the sky background is lower) are also masked, which ex- 
plains the irregular appearance in the comer of the stacked im- 
ages. Fig. [3] shows most of the K s image, resampled 2x2. The 
final image is completely free of any large-scale gradients, and 
the background is perfectly flat except near the brightest objects 
in the field. 

In this release, five stacked images and their corresponding 
weight maps are made available for Y,J,H,K S and YVB118 data 
taken during the first year of public survey operations of the 
Ultra VISTA survey. These images have a zero point of 30.0 AB 
magnitudes for an effective exposure time of one second and 



9 |http : / /casu . ast . ca m . ac . uk/ surveys - pro j ects /vista/ 
technical/astrometric-properties 
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Fig. 3. Full K s mosaic, displayed using a logarithmic stretch. The background level is extremely flat, and is not perturbed near almost 
all bright stars. Several clusters are visible, corresponding to the many rich structures which are present in the COSMOS field. 



a pixel scale of 0.15"/pixel. The weight-maps correspond to 
swarp's image type MAP .WEIGHT which correspond to maps of 
relative inverse variance. Fig.|4]shows an RGB image composed 
K S JY images of a small section of the final field, illustrating the 
excellent image quality and depth of our final stacks. The bright 
saturation limit for stellar sources in these catalogues is ~ 14 
magnitudes in Y and 15 magnitudes in YJHK S bands. 

The images all have a common tangent point, in decimal 
RA, DEC of (1.501163213,2.200973097), corresponding to the 
tangent point of the publicly available IRSA/COSMOS images. 
Each image (uncompressed) is ~ 9Gb in size. This common tan- 
gent point and pixel scale means that the Ultra VISTA survey 
images are pixel-matched to publicly available COSMOS data. 

Finally, to prepare these data products for ingestion in ESO's 
"phase three" system, all the image and table headers produced 
were edited to comply with the Phase 3 requirement document, 
including most of the information which is presented here in the 
FITS header keywords. 

Table [2] summarises the principal properties of each coad- 
ded stack. In each case we report the average seeing over the full 
mosaic, the 95% completeness limit, and the limiting magnitude. 



We also list the typical exposure time per pixel for each stack as 
well as the total on-sky integration time. Summed over all filters, 
this is 55 and 155hrs respectively for the data presented here. 
Seeing on the final stack is characterised using the PSFex tool. 
The average seeing is calculated from a fit to a Moffat| ( 1969) 
profile. We note that in Y band the PSF has slightly broader 
wings compared to redder bandpasses (with a best-fitting Moffat 
/3 parameter which varies from ~ 2.4 in Y to ~ 3.5 in K s ). 



Limiting magnitudes are computed as follows: first, 
SEXtractor is run on each stack using the same detection 
threshold parameters as used for catalogue generation. All pix- 
els belonging to objects to this detection limit are flagged. Next, 
we measure fluxes in apertures of diameter 2" over the entire 
mosaic; any aperture which contains object pixels is discarded. 
The limiting magnitude is then simply computed from the stan- 
dard deviation of fluxes measured in these apertures. Our com- 
pleteness statistics are computed by adding artificial stars to the 
images with average image FWHM and then measuring the frac- 
tion which are successfully detected with SExtractor using the 
same measurement parameters used for the catalogues. 



6 



H J. McCracken et al.: first Ultra VISTA data release 



Table 2. Characteristics of the stacked images. 



Filter 


Typical exposure time per pixel 


Total exposure time 


5cr(2") (±0.1 mag) 


95% comp. (±0.1 mag) 


seeing (")(±0.1") 


Y 


42360 


127080 


24.6 


24.2 


0.82 


J 


49720 


149160 


24.4 


24.2 


0.79 


H 


42520 


127560 


23.9 


24.1 


0.76 


K, 


39400 


118200 


23.7 


23.8 


0.75 


NB118 


23773 


35660 


22.9 ± 0.2 


22.6 


0.75 



Notes. The seeing is computed from a fit to a Moffat ( 1969 1 profile. 
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Fig. 4. RGB image composed of K s , J and Y data respectively. 
The size of this image represents less than l/500th of the total 
area of the field. Sources as faint as K s ~ 22 are easily visible. 



Fig.[5]shows the weight-map from the first year of K s obser- 
vations described in this paper. The intensity at each pixel has 
been converted to an approximate limiting magnitude for a de- 
tection in a 5cr, 2" aperture. It is important to note that our weight 
map is quite uniform, thanks to our adopted observing strategy. 

From these stacks, two sets of catalogues are provided at the 
ESO archive: those extracted on individual images, and matched 
catalogues which use the K s band image as a detection image. 
Aperture magnitudes reported in the catalogues are measured in 
2" and 7.1" diameters respectively. Based on the average stellar 
profiles each of the four broad-band filters, these aperture magni- 
tudes can be "corrected" to pseudo-total magnitudes by adding 

0.35,-0.3,-0.2,-0.2 magnitudes to Y,J,H,K S 2" aperture 

magnitudes. These corrections are not applied to the catalogues 
delivered to the ESO archive but they are applied to the colour- 



colour plots shown in Section 3.5 



3. Data quality assessment 

3.1. Galaxy number counts 

Fig.|6]shows the TTs-band number counts extracted from our cat- 
alogues in comparison with recent literature measurements, in 
particular from the wide-area survey "WIRDS" carried out using 
WirCAM at the CFHT ( |Bielby et al.|20TT) and from COSMOS 
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Fig. 5. Weight-map for the first-year A" s -band data. The intensity 
at each pixel has been converted to an approximate 5cr limiting 
magnitude for an aperture of 2" diameter. The strips at the top 
and bottom of the image have half the average exposure time per 
pixel. 
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Fig. 6. .^-selected galaxy number counts for UltraVISTA, in 
addition to some recent wide-field near-infrared surveys. The 
agreement with previous studies is excellent. 



(McCracken et al. 2010| l. Not surprisingly, our counts agree 
well with the existing COSMOS K s counts but also reach 1 ~ 
mag deeper. We are in good agreement with the other, indepen- 
dent studies covering smaller areas than our work, for example 
dQuadri et al.|2007) . 
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Fig. 7. Difference in position, in arcseconds, with respect to stars 
in the 2MASS, as a function of right ascension and declination 
(upper and lower panels respectively); every second point is plot- 
ted. The solid line shows a running median. The RMS in both 
axes is ~ 0.15". 



Fig. 8. Difference in position in arcseconds for stars between 
the public ACS catalogue of Leauthaud et al. ( 2007) 1 and the 
Ultra VISTA K s stack; every second point is plotted. The solid 
line shows a running median. For both axes the median residuals 
are :£ 0.05". 



3.2. Astrometric comparisons with external catalogues 

We compare the positions in right ascension and declination of 
point sources in 2MASS with those in our Ultra VISTA K s cata- 
logue. This is shown in Fig. [7] Note that, unlike for our photo- 
metric solutions, we do not use 2MASS as our astrometric ref- 
erence catalogue but use instead a densely-sampled catalogue 
from the COSMOS CFHT /-band observations. The absolute 
astrometric calibration of COSMOS is derived from VLA 20cm 
observations (Schinnerer et al. 2004), and these positions are 



known to be offset slightly with respect to 2MASS ( [Capak et al. 
|2007| l, which is indeed what we observe. Our median offsets 
and lcr RMS with respect to 2MASS is (0.00,0.14) arcsec and 
(-0.07, 0.15) arcsec in RA and DEC respectively. 

To verify that our astrometric reference frame is consistent 
with COSMOS, we carried out a similar comparison with stars 
in the COSMOS ACS catalogue ( |Leauthaud eTalJ|2"007l ); this 
is shown in Fig. [8] In RA and DEC, no offset is observed. The 
lcr RMS in both directions for stars selected with 17.0 < K s < 
19.5 is ~ 0.08 arcsec. The internal astrometric accuracy between 
different UltraVISTA bands is expected to be of this order or 
better, i.e., much better than one 0.15" pixel. 

3.3. Photometric comparisons with external catalogues 

We compare the total magnitudes of stars in our catalogue 
(mag_auto) with those in the 2MASS all-sky point source cat- 
alogue ( |Skrutskie et al.||2006| . (Note also that 2MASS is used 
for the photometric calibration of the survey by CASU.) Of 
course, a significant limitation of this comparison is that the 
magnitude range over which sources in UltraVISTA and 2MASS 
overlap is relatively small. Nevertheless, the result of this test is 
shown in Fig. [9] where we plot UltraVISTA-2MASS magnitudes 
for all non-saturated stellar sources and for a total photomet- 
ric error in (2MASS and UltraVISTA, summed in quadrature) 
of less than 0.2 magnitudes. The thick solid line shows a run- 
ning median which is always within 0.05 magnitudes of zero for 
15.0 < mag < 17.0. There is a slight systematic offset visible in 
H (~ 0.03) magnitudes; this could be due to incorrectly rescal- 
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Fig. 9. Difference between total J, H, and K s magnitudes of stars 
in UltraVISTA with sources in 2MASS. The green line corre- 
sponds to a running median. 



ing our exposures to slightly non-photometric images or a real 
offset between the two different photometric systems. 

This Section presents photometric comparisons between 
UltraVISTA and COSMOS JHK S measurements. A large 
amount of near-infrared observations have already been accu- 
mulated on the UltraVISTA field by the COSMOS team. These 
consist of K s - (Mc Cracken et al.1|2010| l and H- band observa- 
tions made with WIRCam on the CFHT and /-band observations 
made with WFCAM on UKIRT. In all cases, these observations 
are shallower than the first-year UltraVISTA data set presented 
here. Since our stacks have the same pixel scale and tangent 
point as the public COSMOS data, to make our comparisons 
we can simply run sextractor in "dual-image" mode, choos- 
ing as detection image the UltraVISTA K s image and as mea- 
surement images the publicly-available COSMOS K s , H and J 
stacks. This approach ensures that no source matching errors are 
introduced. The results of this comparison is shown in Fig. [10] 



8 



H J. McCracken et al.: first Ultra VISTA data release 



RA 

149.4 149.6 149.8 150.0 150.2 150.4 150.6 
1 




0.1 
0.0 

0.1 
0.0 



















• J(RA)|( 



•j • •. • t . ■ . •. v't-- — r~r- 




H (RA) 




0.1 
0.0 



■ J-.UJ JniJ 



J (DEC) 




Fig. 10. Difference in 7, H and total magnitudes between 
Bz,K-selected stars with 17 < K s < 19 in Ultra VISTA and 
COSMOS as a function of right ascension (top three panels) and 
declination (bottom three panels). For clarity, only every fourth 
point is plotted. The thick green line corresponds to a sliding 
median calculated from a window of 100 points. In all cases, 
the differences with the COSMOS photometry is less than 0.1 
magnitudes. 



For test sources we choose Bz/f-selected stars, as described in 
the following section. 

0.1 



In Figs. 10 there is an offset of ~ 0.1 - 0.05 magnitudes 
between Ultra VISTA and the publicly-available COSMOS HK S 
data. We note that at brighter magnitudes, UltraVISTA magni- 
tudes are in good agreement with 2MASS, at least for K s and 
J, bands, and for H the offset reported with respect to 2MASS 
is smaller than the offset with respect to COSMOS magni- 
tudes. There is also some evidence in the K s data of a position- 
dependent offset. Without a third, equally deep data set, it is hard 
to know with certainty the origin of these offsets (especially as 
the VISTA and COSMOS data photometric systems are not iden- 
tical). Furthermore, examining the magnitude of the offsets with 
respect to the COSMOS and UltraVISTA weight-maps they do 
not seem to be correlated with position on the focal planes of 
either instrument (which might be the case if there was a prob- 
lem with the photometric calibration on a chip-by-chip basis). A 
definitive resolution to this issue awaits more involved tests, such 
as photometric redshift comparisons with spectroscopic data, 
which will be the subject of a future article. 




Fig. 11. Seeing FWHM for stars (corresponding to 
SEXtractor's FWHM.WORLD parameter) selected in the BzK 
diagram, as a function of RA and DEC, in the K s stack (every 
fourth point is plotted). As before, the solid green line corre- 
sponds to a running median. The seeing variations are small, of 
order ~ 0.05", and vary principally as a function of RA. 



3.4. Seeing variation across the mosaics 

As described above, the final UltraVISTA stack is comprised of 
six separate "pawprints". At each pawprint the telescope jitter 
displacement is less than the separation between the detectors, 
so no detectors overlap. In general, each OB typically contains 
only images jittered around a single pawprint position, and con- 
sequently the observing conditions, in particular the average see- 
ing is not always identical pawprint-to-pawprint. In first-year 
data presented here, OBs had a mix of maximum seeing con- 
straint between 0.8" and 1.0"; furthermore there is no minimum 
seeing cut. A consequence of this is that when the observations 
are separated paw-by-paw, in some filters, there is a variation of 
around 5% - 10% in average seeing over all 16 detectors from 
paw-to-paw. In the final stack, which is the combination of all 
pawprints, this is visible as bands of regions of slightly different 
seeing. 

In a future UltraVISTA release we will make available stacks 
for which we have carried out a paw-level homogenisation (in 
which each of the six pawprints are convolved by a Gaussian to 
bring them to a common FWHM). For the moment, we report 
here that this effect is important for the K s and H- band stacks. 
In Fig.[TT]we show the seeing, calculated from SExtractor's 
FWHM_WORLD parameter (which is derived from the isophotal area 
of the object at half maximum, and so may not be comparable 
to the figures listed in Table [2}, as a function of right ascension 
and declination. Because of a sequence of pawprints with signif- 
icantly better seeing, there is around a 5% variation in seeing as 
a function of right ascension. 

3.5. Colour-magnitude and colour-colour diagrams 

The large number of sources in our catalogues combined with 
our excellent seeing and high signal-to-noise means that we can 
investigate in detail the distribution of objects in colour-colour 
space. In Fig. 12 we plot the (J-K s ) vs K s distribution of sources 
in our K s selected catalogue. The stellar locus is clearly visible as 
a narrow ridge of constant (J-K s ) colour (which one can confirm 
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Fig. 12. Two-dimensional histogram showing ( J — K s ) corrected 
aperture colour as a function of K s total magnitude; the grey level 
at each bin in magnitude-colour space corresponds to the surface 
density of objects. The narrow ridge clearly visible at {J - K s ) ~ 
-0.2 corresponds to the location of stellar sources. 
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Fig. 14. Stellar locus for bright and faint stars in Ultra VISTA 
(shown as points and dots respectively) in the (B - z) vs (z - 
K) corrected aperture colour-colour plane. Bright and faint stars 
occupy the same location in colour-colour space. 




Fig. 13. Two-dimensional histogram of (B - z) vs (z - K s ) cor- 
rected aperture colour for UltraVISTA. All sources detected to a 
5cr limit in K s auto magnitudes are shown. The stellar locus is 
clearly visible as a ridge at blue (z - K s ) colour. 



I 




-0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 

(y-j) 

Fig. 15. Two-dimensional corrected aperture colour-colour (H - 
K s ) vs (Y - J) histogram for all sources with total magnitude 
16.0 < K s < 23. Stars and galaxies are cleanly separated. The 
stellar locus corresponds to the 'bump' visible at 0.2, 0.1 in (Y - 
J) vs (H - K s ). 



by overplotting on this diagram the location of stars identified in 
the ACS catalogue). 

Next, we consider the distribution of objects in optical and 
near-infrared colour-colour space, turning first to the "BzK" di- 
agram as this allows us to cleanly separate stars and galaxies. 
We use the publicly-available COSMOS B and z Subaru images 
( |Capak et al.| [2007) and transform the B and z magnitudes in 
each catalogue following the recipes in McCracken et al. ( 2010|> 



almost parallel to the stellar locus but redder in (B - z) is a sec- 
ond "ridge" which is comprised mainly of lower-redshift passive 
galaxies ( |Lane et al.|[2007) |Bielby et al.1[20TT] l. Thirdly, the di- 
vision between lower-redshift normal and star-forming galaxies 
(the "sBzK" galaxies of Daddi et al. ( 2004| l) is clear. 



to bring our system to the "BzK" system defined in|Daddi et al. 
( |2004] >. 

The result is shown in Fig. 



13 as a two-dimensional grey- 



scale histogram; in this diagram and all subsequent diagrams 
we show all objects detected to 5<x in ^ s band aperture magni- 
tude. Several interesting features are clearly visible in this di- 
agram: firstly the stellar locus, which is apparent as the long 
"ridge" feature which is relatively blue in (z - K s ); secondly, 



In Fig. 14 we show a magnified view of the stellar locus in 
the BzK diagram, and we show both bright and faint stars. The 
position of the stellar locus does not depend on the magnitude 
limit, which demonstrates that there are no magnitude-dependent 
effects present in our data which could arise if there were issues 
related to an incorrect sky-subtraction. 

We also consider the distribution of galaxies in the purely 
near-infrared colour-colour space (H - K s ) vs (Y - J), shown in 
Fig. [15] Again, stars and galaxies are cleanly separated. 
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4. Summary and conclusions 

In this paper we have described the first public UltraVISTA data 
release. This data set comprises five high-quality image stacks 
representing a unique combination of depth and areal coverage at 
near-infrared wavelengths. Our stacked images reach 5<x depths 
in aperture of 2" diameter of ~ 25 in Y and ~ 24 in JHK S bands. 
Furthermore, it is worth noting that these depths are in agree- 
ment with the expected sensitivity of the telescope at the time 
of writing the original UltraVISTA survey proposal. To these 
limits, our K s catalogue contains 216,268 sources. The lcr as- 
trometric RMS in right ascension and declination for stars se- 
lected with 17.0 < K s < 19.5 is ~ 0.08 arcsec in comparison 
to the publicly-available COSMOS ACS catalogues. Each of the 
stacks has sub-arcsecond seeing and the FWHM variation over 
the images is less than 5% in most bands. Our number counts 
and photometric calibration are in good agreement with previ- 
ous studies. 

The images and catalogues described here are publicly avail- 
able from the ESO archival 

At the present time of writing (April 2012), a further 250 
hours of UltraVISTA observations have been completed. We in- 
tend to deliver regular releases of Ultra Vista data products as 
the observations proceed towards the total 1800h of observation 
time allocated to the project. 
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